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Abstract 

QQ ■ We investigate the reaction pn — )■ dinT in the framework of Chiral Perturbation Theory. For 

the first time a complete calculation of the leading order contributions is presented. We identify 

various diagrams that are of equal importance as compared to those recognized in earlier works. 

^ ' The diagrams at leading order behave as expected by the power counting. Also for the first time 

the nucleon-nucleon interaction in the initial, intermediate and final state is included consistently 

►^ , and found to be very important. This study provides a theoretical basis for a controlled evaluation 

QQ I of the non-resonant contributions in two-pion production reactions in nucleon-nucleon collisions 

m ' 

o 
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I. INTRODUCTION 

As the first strong inelastic threshold of the nucleon-nucleon system, pion production in 
A^A^ collisions has attracted a large number of theoretical as well as experimental works. 
Due to advances in the experimental methods measurements in the threshold region became 
possible in the beginning of the 1990ies, and it quickly became clear that the theoretical 
models existing at the time [l| were not able to describe these data — see Ref. [2] and 
references therein. They fell short by a factor of two for reactions with an isoscalar A^A^ pair 
in the final state, while there was a discrepancy of even more than an order of magnitude 
for those with an isovector A^A^ pair in the final state. The numerous efforts to improve the 
phenomenological approaches 3|, \^ , although quite successful for various observables |5|] , did 
lack for systematics, gave contradictory answers, and involved a sensitivity to unobservable 
short range effects. 

To overcome those deficiencies, in recent years considerable theoretical efforts went into 
the development of an effective field theory that can be applied to the reactions A^A^ — )• NNtt. 
In early studies it became clear, however, that the original powercounting by Weinberg [a, |7[ 



9j (see also Ref. 
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where 



needs modifications in order to arrive at a convergent expansion 
it was pointed out that the naive power counting using the heavy baryon formalism is not 
applicable above the pion production threshold — the necessary adaptions are outlined in 
Ref. [ll|). Indeed, for neutral pion production in pp collisions, the corrections due to the 
next-to-leading order (NLO) increased the discrepancy with the data and, moreover, the 
next-to-next-to-leading order (NNLO) contributions turned out to be even larger than the 

n 

NLO terms [12] . The origin of these difficulties was identified quite early by Cohen et 
al. [l3|, see also [l4|, who stressed that the additional new scale, inherent in reactions of 
the type of A^A^ — )■ NNn, needs to be accounted for in the power counting. Since the two 
nucleons in the initial state need to have sufficiently high kinetic energy to produce the 
on-shell pion in the final state, the initial center-of-mass momentum needs to be larger than 

0.4, (1) 

where ttItt and Mn refer to the pion and nucleon mass, respectively, and A,^ denotes the 
typical hadronic scale, here identified with M^ to get a numerical estimate for the expansion 
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parameter. The proper way to include this scale was presented in Ref. 15| and implemented 



in Ref. 



Q 



16| . see Ref. [2j for a review article. As a result p-wave pion production is governed 



by tree-level diagrams up to NNLO 15|, [iTJ] in the modified power counting scheme. On the 
other hand, for s-wave pion production, pion loops, studied in detail in Ref. 18|, start to 
contribute already at NLO. The individual loop contributions turned out to grow linearly 
with increasing A^A^ relative momentum which resulted in a large sensitivity to the employed 
A^A^ wave functions 19| . However, it was demonstrated in Refs. 18|, |20(] that all irreducible 



loop contributions at NLO cancel altogether which was needed to keep the scheme self- 
consistent. Furthermore, the proper separation of irreducible terms at NLO from reducible 
ones in the loops resulted in an increase of the most important LO operator for charged pion 
production, first investigated in Ref. pj, by a factor of 4/3. This increase was sufficient to 
overcome the apparent discrepancy with the data in the reaction pp — )• dn^. The neutral pion 



channel is more challenging — it still calls for a calcu^ 
First steps in this direction were taken in Refs. 11 



ation of subleading loop contributions. 



211. We further emphasize that the 
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A(1232) isobar should be taken into account explicitly as a dynamical degree of freedom [13| 
because the Delta-nucleon mass difference, AM, is also of the order of Ptjn- This general 



m 



23|, see also 



argument was confirmed numerically in phenomenological calculations 
Refs. |l6l . |20| where the effect of the A in NNtc was studied within ChPT. 

The progress in the theory of isospin conserving pion production in nucleon-nucleon 
collisions allowed us to perform a complete leading-order calculation of charge symmetry 
breaking effects in pn — > dn^ 2J|, see also Refs. 25|, |26| for related works. 

Given the developments just discussed, one now should be in the position to investigate 
two-pion production in A^A^ collisions using the same formalism. Such calculations do 
not exist yet, although in the pioneering works of Refs. 27|, |28|] some constraints from 
chiral symmetry were already implemented. A more recent phenomenological analysis of 
this reaction can also be found in Ref. 29j. For two-pion production in A^A^ collisions at 



threshold, the center-of-mass momentum of the initial nucleons is necessarily larger by a 
factor of \/2 compared to the single-pion production. Thus we now have 
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2m„MN ~ 510 MeV . 
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With this large momentum scale, the value of the expansion parameter is 

0.54 
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and one may question the applicability of CliPT for the two-pion production case. However, 
while the expansion may converge slowly, it is still meaningful at least close to the chiral 
limit and therefore should also provide some guideline on the hierarchy of diagrams. Thus 
it is still interesting and important to investigate the structure of diagrams that contribute 
at the lowest orders. Note that, since the NN momentum at the threshold for two-pion 
production is quite close to the one for single-pion production, for the purpose of power 
counting below we identify the two. 

The amount of experimental data for A^A^ — )■ NNtttt has increased greatly in recent 
years 30l-l33| and even a polarized measurement exists 3J]. However, up to now there are 
no data available sufficiently close to threshold to allow for a straightforward comparison 
with our leading order calculation of the reaction pn — )■ dnn. To extend our calculation 
to higher energies the inclusion of the A as an explicit degree of freedom will be neces- 
sary (it enters at NLO). In addition, at NNLO, besides a large number of loops, also the 
first A^A^ — )■ NNtttt counterterms start to contribute. The physics of those is dominated 
by heavier baryon resonances — most prominently the Roper resonaiice, as follows from 
previous phenomenological studies of the reactions A^A^ — )■ NNtttt 27H29| . Two pion pro- 
duction and the role of the Roper resonance, in particular, has been studied extensively 
in the single-nucleon sector, namely in the reactions 7A^ — > ttttN and vrA^ — t- ttttN both 



phenomenologically, see, e.g., Refs. 



ChPT studies |39|, 



35 



and within ChPT |39|-|43| q In particular. 



B 



40| predict that the double-neutral-pion photoproduction is significantly 



enhanced near threshold as compared to the production of the charged pions due to the 
contribution of chiral (pion) loops at NLO (0(p^)). This prediction was confirmed experi- 
mentally in Ref. 4J]. Also it was found in Ref. |39| that the A-isobar contribution, which is 
potentially important because the A— nucleon mass difference is of comparable magnitude 
to 2r7i^, turns out to be insignificant due to specific cancellations of the individual contri- 
butions. The Roper resonance contributes to this process at one order higher, i.e. at N^LO 
(0(p^)). The effect of the Roper resonance encoded in low-energy constants was found to 
)e sizable although significantly smaller than the contribution of the leading chiral loops 
40| . Similar conclusions about the role of the Roper resonance were drawn based on studies 
of ttN —7- ttttN. In this case, however, the loops were found to be negligible 42| and the 



^ Note that for the reaction irN — >■ ttttN within ChPT we refer to the calculations at order 0{p^) only, see 
references cited in 39l-l43l| for numerous other studies. 



leading contribution was provided by tree level terms. For two-neutral-pion production in 
A^A^ collisions we expect the production mechanism to be in accord with the dimensional 
analysis, i.e it should be dominated by the tree level graphs at LO. However, an explicit 
study of loops at NLO will be necessary to support our uncertainty estimates. 

The paper is structured in the following way: In Sect. II our theoretical framework is 
introduced. Results for the reaction pn — )■ (Itt^tt^ are presented and discussed in Sect. III. 
A summary and some concluding remarks are provided in Sect. IV. Technical details of our 
calculations are summarized in two Appendices. 

II. THEORETICAL FORMALISM 



Our calculations are based on the effective chiral Lagrangian 45N47l|: 



+^N^T ■ a ■ f Vtt + 7r^7r(7r ■ Vtt) ) A^ + ■ ■ ■ 
£(1) = ^ (iN^T ■ (tt X Vtt) ■ VN + h.c.) - ^ f/ ^ [iNW ■'ka-VN + h.c] + 

'^2 + C3 - ^^ TT^ - C3(V7r)^ - 2cim^7r^ - -(c4 + —T-r-yij^.eahcCyk'TcdiT^adjT^h 



P 



N 
+ ■■■ (4) 



where the superscripts and 1 denote the leading and next-to-leading order Lagrangian, 
respectively, f-,^ denotes the pion decay constant in the chiral limit, qa is the axial- vector 
coupling of the nucleon, A^ and vr are the nucleon and pion field, respectively, and a (r) 
denotes the Pauli matrix in spin (isospin) space. The ellipses in the above equations represent 
terms which are not relevant for the present study. 

As mentioned in the Introduction, the system is characterized by three scales, where we 
use 

m^ < p ~ plhl ~ p£i ~ AM < A^ . 



The counting rules here are similar to those for one-pion production used in 15l-ll8| (for a 



review see also 2j). For the Weinberg- Tomozawa (WT) ttN — ?■ ttN vertex, as discussed in 
detail in Ref . [l8| , the leading term should be proportional to 2ut, where Uj, is the energy of 
the outgoing (on-shell) pion. The diagrams that contribute at the lowest order are shown in 
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FIG. 1. Leading order diagrams for the reaction pn — )• dTrvr at threshold. Sohd (dashed) hnes 
denote nucleons (pions), filled ellipses correspond to continuum NN wave functions (including 
plane wave) in the initial and intermediate state, the outgoing black line denotes the deuteron. 
Subleading vertices are marked as 0. 



Fig. [U and start to contribute at order x ~ itIt^/M^. Obviously, the number of diagrams is 
significantly larger for two-pion production than for one-pion production already at LO — 
after all the latter appears as building block of the former in diagrams a)-d). For Fig. [T]i), 
we note that the naive dimensional analysis shows that this diagram contributes at next- 
to-leading order. However, due to the existence of the two-nucleon cut in the intermediate 
state, this diagram is promoted to leading order following the very same logic adopted by 
Weinberg [3] . In addition. Fig. [T] e)-i) show genuine two-pion production diagrams. Note 



that the only non-resonant diagrams included in the earlier studies of Refs. [27|, |28| are those 
depicted in the second line of Fig. [H whereas in Ref. 29[ such diagrams are not considered 
at all. 

Another potentially important class of contributions arises from diagrams where one of 
the nucleons is excited to the A(1232)-resonance. However, in the reaction considered here 
the A^A system can only couple to the A^A^ system in the intermediate state. The A^A^ 
systems in the initial and final states have to have isospin zero, cf. below, while the A"A 
system has always isospin one. The A(1232) can only appear when the first pion is emitted 



and, thus, the power counting for the A(1232) in pn — )■ (ivrvr is similar to the one in single-pion 
production. In particular, the A(1232) starts to contribute at NLO due to the fact that the 
A propagator is suppressed by 1/p as compared to l/mT^ in the nucleon case. Moreover, the 
aforementioned isospin arguments prevent the potentially dangerous situation when there is 
an extremely enhanced A-propagator (AM — 2r7i^ — recoil)"^ in the A^A intermediate state 
that would occur from a A excitation prior to the first pion emission. It is interesting to 
note that a similar observation was made in Ref. J39| in the reaction 7A^ — t- nnN where the 
corresponding term (AM — 2771,^)"^ did not show up in the denominator because of an exact 
cancellation with the same term in the numerator. 

Once the production operators are constructed, they should be convoluted with A^A^ wave 
functions to account for the nonperturbative character of A^A^ interaction [7]. In principle, 
the A^A^ wave functions should be calculated on the same basis as the production operators. 
However, up to now ChPT A^A^ interactions are only available for energies below the one- 



pion production threshold 



48 



g|-l50| . Therefore, in this work we adopt the so-called hybrid 
approach that consists in the convolution of the production operator calculated within ChPT 
with phenomenological A^A^ wave functions from realistic A^A^ models, in particular from the 
CCF model |511]. We utilize also wave functions generated from the CD-Bonn A^A^ potential 



52| | so that we can examine in how far our results depend on the specific choice of the A^A^ 



wave functions. The A(1232) resonance is included as explicit degree of freedom in the 



construction of the CCF potential [5l|] as well as in an extended version of the CD-Bonn 
potential 53|. This will allow us to extend in a future work the present calculation to higher 
orders in a straightforward manner. The ellipses in the Feynman diagrams of Fig. [T] refer 
to the A^A^ wave functions and represent either the Initial State Interaction (ISI) or the 
Intermediate State Interaction (ImSI) of nucleon pairs supplemented with the plane wave 
terms. 

Based on the Lagrangian given above, it is straightforward to calculate the amplitudes 
corresponding to the diagrams shown in Fig. [H Furthermore, since we are only interested in 
the energy region very near to threshold, only a few partial waves contribute to the reaction 
amplitude. At threshold, the vrvr system is in a relative S-wave and its isospin can be either 
zero or two as required by the symmetry of a system of identical bosons. However, since the 
isospin of the initial A^A^ system can only be zero or one and the deuteron is an isosinglet 
the isospin of the produced nn pair has to be zero and, consequently, also the isospin of the 



TABLE I. The values of the amphtudes for the individual diagrams (in units of 10 ^xMeV ^] 
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initial A^A^ state is zero. Then the conservation of parity and angular momentum further 
requires that the total angular momentum of the system should be one and the NN pair 
in the initial state can only be in the ^Si or ^Di partial waves. The Pauli principle is also 
satisfied in this case. The explicit expressions for the amplitudes of these two partial waves 
can be found in Appendix A. 

III. RESULTS AND DISCUSSION 



As discussed in the introduction, some leading order diagrams depicted in Figs. [T^)- [T^) 
are missing in Ref. |27|. First, we want to discuss the role of the diagrams [It)- [Hi). As 



TABLE II. Values of the individual amplitudes for the Si partial wave (in 10 xMeV ). 
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can be seen from the figure, these diagrams are closely related to the one-pion production 
case. There the leading order contributions come from the direct pion emission and the 
rescattering term. Thus, one may expect that the amplitudes [I^)- [Hi) reveal a pattern 
similar to the one-pion production case. For example, one may guess that the amplitude of 
Fig. [1^) should be much smaller than the one of Fig. [T]d) due to the fact that for single-pion 
production the rescattering term dominates over the direct one by an order of magnitude. 
However, it should be noted that in our case the NN amplitude in the intermediate state is 
fully off-shell and it turns out that the real part of diagram b) is larger only by a factor of 
about 3 as compared to the one of diagram a) — cf. Table [B In addition, due to the fact 
that the A^A^ amplitudes for the ISI appear in coupled channels (^S\-'^D\) the production 
amplitudes, as given in Tabled do not have the same phase. This fact makes the interference 
between the individual amplitudes very different from the one-pion production case and a 
direct comparison becomes difficult. The role of the NN interaction in the initial and 
intermediate states can be understood from looking at the individual contributions to the 



TABLE III. Values of the individual amplitudes for the Z?i partial wave (in 10 xMeV ). 
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two-pion production amplitudes listed in Tables [TTland lllli We find that the values of these 
contributions are basically comparable with each other, as expected from the power counting. 
The values of the amplitudes of Fig. [T]d) and Fig. [TJi) are about a factor of two to three 
larger than the other two amplitudes [1^) andlH^), which is still a reflection of the dynamics 
governing the one-pion production. We observe that the resulting values of the aniplitudes 



[T^)- [Hi) are similar in size to those of [Tf)- [!}) considered in previous studies (27|, l28| and, 
thus, are important. We also find that the A^A^ interaction in the initial and intermediate 
states plays an important role in the calculation, especially for the diagrams [T^)-[T]i). In 
particular, one can find from Tables HTl and IIIII that the inclusion of the A^A^ interaction in the 
ISI or the ImSI generally reduces the magnitude of the amplitudes as compared to the Born 
amplitude and that the cancellation among the individual amplitudes of each particular 
diagram is significant. The largest cancellation takes place between the amplitudes with 
and without the ImSI, i.e. between (^q and Aimsi) and {Aisi and Aisi-^imsi) where 
the amplitudes are defined in Appendix A, see also Tables [III and [1111 ■ For example, the 
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value of the S-wave amplitude of Fig. [T^) becomes a factor of two smaller than the Born 
amplitude after including the ISI. When the ImSI is also considered, the full amplitude is 
further reduced so that it is finally a factor of six smaller than the Born amplitude. This 
observation shows explicitly that the ISI and the ImSI are very important quantitatively. 
Since they also strongly influence the phase of the amplitudes, their proper inclusion is 
compulsory, especially for studies of polarization observables. 

The strength of the contact terms, that appear in Fig. [T^), is given by the low energy 
constants Cj that can be determined from vrA^ scattering. Because we are only interested 
in the amplitudes at threshold, we only need the values of ci, C2 and C3. In this work, we 
take the values from Ref. 5J| where the values of q's are obtained by fitting vrA^ threshold 
parameters. To the order we are working, three solutions are offered in that paper. One 
corresponds to the results without considering the A (1232) explicitly, and the other two 
are results where the A(1232) is included but with different choices of the tcNA coupling 
constant /i^. It is well known that C2 and C3 can be largely accounted for by the contribution 
from the A(1232) 55|, and therefore the values of C2 and C3 can change significantly between 
extractions with and without inclusion of the A(1232). One may expect that using different 
sets of parameters will affect the value of the amplitudes of Fig. [1^) strongly. However, it is 
interesting to note that the amplitude of Fig. [1^) is independent of those choices, because 
it only depends on the value of the linear combination 

2mlci + g?g°(c2 + C3 - (^^/(SMjv)) , (5) 

with q^ = q2 = ni.^ for the pion energies in the kinematics relevant here. From Table 1 of 
Ref. |i5j] one can see that the values of ci and C2 + C3 are not affected by the A(1232) up 
to order 0{p^). Thus, although the values of C2 and C3 can significantly change individually 
by considering the A(1232) explicitly, the sum of them is not affected up to order O(p^), 



as shown explicitly in Ref. [55!]. Hence at LO the contribution from contact terms are not 
influenced by the A(1232), which is also consistent with the power counting used. It is also 
interesting to note that the linear combination of low energy parameters displayed in Eq. ([5]) 
is large, for the individual terms interfere constructively — this is in contradistinction to 
ttN scattering, where q^ = — g^ = ^^n holds at threshold and then physics is governed by 



the very small isoscalar scattering length [55|, see Ref. [56| for a recent accurate extraction 
of a~^ from pionic atoms. Thus, once the reaction pn — t- dn^n^ can be studied with sufficient 

11 



accuracy, it might be a good source of information for the Cj individually. Note also that 
in Ref. J5J] the l/M^r corrections are not included explicitly in the fits to nN data but 
absorbed effectively into the LECs Cj. To be consistent with their treatment we also drop 
the I/Mtv term in Eq. (|5]). Furthermore, the residual combination C2 + C3 which enters 
Eq. (|5]) is totally determined by the S-wave threshold parameters. Numerically the result 
for diagram [T^) turns out to be one of the largest individual contributions — cf. Table HI 

The diagramcl of Figs. Hf)- [T}) are included in the previous study of Ref. j271] , except for 
the initial state interaction which is not considered in that work. However, as can be seen 
from Tables HTl and IIIH the consistent inclusion of the ISI changes the value of the amplitudes 
[T]F)-[T}) significantly and is thus important. 

Besides the interference between the individual amplitudes of each particular diagram, 
the interference between the contributions from different diagrams is also very significant as 
can be seen from Tables [TTl and IIIII In order to examine the sensitivity of our results and 
especially of the interference pattern to the employed A^A^ interaction we also performed 
calculations with a different A^A^ model, namely with the CD-Bonn model \b2\. We found 
that the values of the individual amplitudes vary in a reasonable range (up to 30%) due to 
differences in the A^A^ wave functions whereas the sum of all amplitudes varies just by a 
small amount (around 10%). This variation should be partly absorbed by the contribution 
of the {N^NfTx-K contact term at N^LO. 

To check the convergence of chiral expansions based on counting rules, it is necessary 
to calculate explicitly subleading contributions. As those calculations are not available yet, 
we may try to compare the results from the leading order contributions to experimental 
observables: our result should be within 50% of the data, given the relatively large value 
of the expansion parameter. The applicability of our s-wave calculation is restricted to the 
near-threshold regime and, therefore, data close to the threshold are needed for a sensible 
comparison of theory with experiment. Unfortunately, at present the empirical information 
for the (ivrTT channel is very scarce. Only data down to excess energies Q = 70 MeV {Q = 
y/s — y^Sthi) are available and the uncertainties are still very large. Thus, it is difficult to 
draw any conclusions about the convergence in the present work. In order to provide some 
comparison of our calculation with the existing data, we assume that the matrix element 



We checked, using the methods of Ref. 



Ill ], that the sum of the contributions from Figs. [T^)- HJ) is 



independent of the choice of the pion field, as required by field-theoretic consistency. 
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Q(MeV) 

FIG. 2. Result of our ChPT calculation for the cross section of pn — )• d-K^-K^ based on the assumption 
of a constant production amplitude. The shaded area corresponds to the uncertainty estimate based 
on a dimensional analysis. 

is constant and therefore take the energy dependence as originating only from the three- 
body phase space. Our result is roughly a factor of two smaller than the central value of 
the data point at Q = 70 MeV. We want to stress, however, that the uncertainties of the 
experiment in question is quite large (even a factor of 25 larger than the data themselves). 
For future reference, in Fig. [2] we show the central value for the cross section predicted 
under the assumption that the matrix element is a constant and that no higher partial 
waves contribute. We also indicate the uncertainty band based on a dimensional analysis. 

IV. SUMMARY 



We presented a complete leading order calculation for the reaction pn — )■ d-K^ir^ at thresh- 
old within chiral perturbation theory. There is no free parameter in the calculation at this or- 
der. We included various additional diagrams as compared to previous investigations |27H28| 
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and we also made several technical improvements. Our most important findings, summarized 
in Table HI are 

• all diagrams evaluated are of similar importance; 

• there are sizable interferences between the individual contributions; 

• the accurate inclusion of the NN interaction in both intermediate as well as initial 
states is very important. 



We also stress that for two-pion production the expansion parameter 

0.54 , 



X 



(2) 
Pthr 



M, 



N 



(2) 

where p^hr denotes the initial NN momentum at the two-pion production threshold, is rather 
large. It is therefore important to calculate higher order contributions to check the rate of 
convergence. We argue in this paper that despite of the proximity of the A(1232) to the vrvrA^ 
threshold the potentially most important A^A intermediate state is not allowed for the reac- 
tion pn — > dnn because of isospin conservation. Intermediate states containing the A(1232) 
must be of the kind vrA^A, i.e. can only occur after one-pion emission. Therefore, the role 
of the A(1232) resonance in the reaction considered is expected to be analogous to that 
in one-pion production. In particular, the A(1232) starts to contribute at next-to-leadin; 



order. Furthermore, it is known from phenomenological studies of A^A^ — )■ NNnn 27|, |28| 
that the Roper resonance can play a significant role already near threshold and that it will 
become even more important when considering a larger range of energies. However, this res- 
onance is not included explicitly in the present study. One may expect that its contribution 
is absorbed into some low energy constants. The contribution of the Roper resonance to 
the Ci parameters, that scale the strength of the leading isoscalar vrN scattering, has been 
discussed, e.g., in Ref. [55|. It seems that it plays only a minor role here. Similar conclu- 
sions were drawn from systematic studies within ChPT of the double-pion photoproduction 
40| and the reaction vrA^ — )■ vrvriV 4ll-l43| near threshold. In particular, for 



process [39 



the reaction 7p — !■ ti^tt^p the contribution of the Roper was found to be rather moderate 
as compared to the large contribution of chiral loops 40|]. A recent model calculation of 
11 N —7- ttttN by S. Schneider et al. [35| suggests also that the Roper resonance plays a rather 
minor role. On the other hand, the Roper might contribute significantly to the (A^^A^)^7r7r 
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counterterms, which enter at NNLO in a ChPT calculation of A^A^ — )■ NNnn. Based on 
the discussions above, one can not expect that our current calculation can describe the 
experimental data well at higher energies. However, the presented calculation provides an 
estimate for the contribution of the non-resonant background near threshold. It therefore 
forms a basis for future studies and is thus a precondition to extract reliable information on 
the Roper resonance from near-threshold experiments. 
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Appendix A: Reaction amplitudes 

In this appendix, we present expressions for the amplitudes that we consider in this work. 
To obtain the partial wave amplitudes, we follow the technics developed in Ref. |57|. 

Eor the reaction A^A^ — )■ (Itttt at threshold, the initial state can only be in the ^Si or 
^Di partial waves and the total isospin is zero. Thus, the amplitudes for A^A^ — t- cItt-k at 
threshold can be written in the general form 

M = {t. S^''^U% + t . S^'^^Ufj,) ■ (vf . vf)Xo , (Al) 

where iS| = X2''^2:5i^i ^^^ '^i ~ K'^*-?' ~ 3^*^i) ' xf^'^aCTjXi denote the normahzed 
spin-orbit structure of the initial nucleons for the ^Si and ^Di partial waves, respectively, 
Xq = (p2^^i denotes the isospin structure of the initial states, Aj^n and Aj^^i are the 
corresponding partial wave amplitudes, vf is the isospin wave function of outgoing pions 
and e'is the deuteron polarization vector. Here Xi(X2) and ipi{ip2) are spinors of the initial 
nucleons in spin and isospin space, respectively, and n denotes the unit vector of the relative 
momentum of the initial nucleons. The expressions for Aj^ii are obtained by projecting 
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the amplitudes that contribute to a given order into corresponding partial waves using the 
technics developed in Ref. [57|. To leading order, the relevant diagrams in our work are 
shown in Fig. [1] Each diagram in Fig. [1] represents a set of diagrams, in which possible 
A^A^ interaction in the initial and intermediate A^A^ states are considered. For example, the 
partial wave amplitudes of Figs. [T^)-d) originate from several parts as follows: 

aPW _ aPW , aPW : aPW : aPW ^ a 9) 

v^ — ..^Q -t- ~^jsi ^ -^ImSI ^ -^ISI+ImSI ' l^^J 

where the index '0' stands for the Born diagram, and 'ISI', 'ImSI' or 'ISI+ImSI' correspond to 
the initial state interaction, intermediate state interaction, or both initial and intermediate 
state interactions included in the amplitudes, respectively. 'PW represents the two possible 
partial waves in the initial state, which can be either ^Si or ^Di. For Figs. [I^)-i) the 
partial wave amplitudes contain only two parts, namely Aq^ and A^^j because there is no 
intermediate state interaction. 

To compute diagrams with the A^A^ interaction in the initial and intermediate states we 
take the A^A^ scattering amplitudes from some potential models. In our calculation, instead 
of the commonly used T matrix, the Ai matrix is used. These quantities are related by 
M = -%Ti'^M%T. For the initial state interaction M^^^^' denotes the A^A^ half-off- 
shell M. matrix, where PW — )■ PW represents the transition from the partial wave "Piy" 
to ''PW " with "PPF" and ''PW " being ^^i or ^Di in our case. For the intermediate 
state interaction, the only possible transition is ^Pi — >^ Pi, so we just use Aij^]^ ^ to 
denote the fully off-shell A^A^ M. matrix. For the deuteron-NN vertex we adopt the notation 

]. We use u{p) and w{p) to denote the A^A^ components 



and structures used in Ref. 



of the deuteron wave function corresponding to S-wave and D-wave respectively. Clearly, 
the bound (deuteron) wave functions u{p) and w{p) should be calculated using the same 
potential model as the A^A^ scattering amplitudes in the continuum state. In this work we 
use the CCF model [511 to generate the A^A^ amplitudes and the deuteron wave functions. 
Since the AA channel is considered in the CCF model, that can also couple to the deuteron, 
the A^A^ part of the deuteron wave functions, u{p) and w(p), are normalized as: 

/ 10^ (""^^^^ + '^^^^') = °"^^^^ • ^^^^ 

In this work, we adopt the following values of the parameters: /^ = 92.4 MeV, qa = 1.32, 
m^ = 139.58 MeV and Mjsi = 938.27 MeV. The loop integration is regularized using the 
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cut-off method with pcut = AchPT ~ 1 GeV. We checked that the dependence of the results 
on this parameter is small - the results change by 10% when A is increased to 10 GeV. 

1. Amplitudes of Fig. [1] a)-d) 




/ [ 



ii 




-.{ 



4/ 
/ 



P -P P -P 

U) (B) 



ii 



P -P 

(C) 



P -P 



FIG. 3. Building blocks for construction the diagrams a)-d) of Fig. [T] 



To get the partial wave amplitudes of Fig. [T] a)-d) it is convenient to separate these 
diagrams into some building blocks at the position where the intermediate state interaction 
may occur, as shown in Fig. [3l Here the diagrams (A) and (B) correspond to the transition 
from the ^Pi partial wave that is realized in the intermediate state in the diagrams of Fig. [1] 
a)-d) to the ^^i deuteron state. Correspondingly, the diagrams (C) and (D) represent the 
transition for the ^5*1 — ^ Di initial state of the two-nucleon pair to the intermediate ^Pi 
state. The full amplitudes can then be obtained from the expressions for these building 
blocks by supplementing propagators, possible NN amplitudes and loop integrations. The 
amplitudes of Fig. Wi,A) and (B) can be expressed in general as: 



M = {e'- 5^'^^] Xi 



TT 



B 



(A4) 



where ttj is the isospin wave function of the vr meson in the final state. Here we introduced 
Sl = 2^vkPkX2'^20'jXi ^-iid Xi j = (p2T2y^(pi to dcuotc the normalized spin and isospin 
structures of the intermediate NN state, and p = p/p is the unit vector of the relative NN 
momentum in the intermediate state. The expressions of B for diagrams (A) and (B) are : 



BAip) 
Bb{v) 



3 — u — n^'^^^^^r 



(A5) 



4 gAm^M^ f Pdl f dflp dfli 



Vs n 



1 
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fw.)-f: 



-^wip)ip.ir+i=iuip) + ^)p.i 



(A6) 



Correspondingly, the expressions of the building blocks for Fig. Et^C) and (D) can be given 
in the following way: 



where B ^^ and B ^^ for the diagrams (C) and (D) read: 



(Xl ■ 7f)Xo 



BT{P) 
Bf{p,l) 



__P_gArrhr_ 

1 QAf^-K f dQp dQi pp ■ I — I 






in 4n (p*_/)2 + ^2 



dVtpdVLil{p-if 



2p 



p ■ I 



(A7) 

(A8) 
(A9) 

(AlO) 
(All) 



4 f^ J An An (p-- /)2 + ^2 

Given the building blocks shown above and the corresponding A^A^ scattering amplitudes, 
it is straightforward to construct the amplitudes of Fig. [T]a)-d). For example, the amplitudes 
corresponding to diagram [T^) are: 



a: 



a,PW 




Ns ■ B^^ 



(Pi 



A 



a,PW 
ISI 



Ns- E 

PW 



Pi 
^2, 



■ BAiPi 



pI 

fdp A;l^^-^^'(p,,p) 

27r2 



(A12) 



M, 



N 



{2M^Y 



pA 



P2 -«e 



nPW 



{P) 



pz _ p| 



te 



Ba{p) 



Aa,PW TIT rsPW / 
ATL<^r = Ns-Bc {Pi 



M, 



N 



A 



^ImSI 



a,PW 
ISI+ImSI 



AT; 



PW 



Pi - Pi 

p^dp M 



l^dl M'^V'^'^iPiJ) 



M, 



N 



27r2 



(2M, 



N 



P 



PW^PW 

NN 



'{Pi,P) 



M, 



N 



27r2 



{2Mn 



IHI A^^V'^^Hp,/) 



2 



le 



P P2 ~ "^^ 
Sa{1), 



p{ — le 



BE'^'ip) 



(A13) 
Ba{1) , (AM) 



M 



N 



pz 



Pi — le 
(A15) 



27r2 {2MnY P - pI 

where p2 = p[l^,. and pi = pl^^. Here A^^ is the symmetry factor which is obtained by 
considering the interchange of identical particles in the initial, intermediate and final states. 
This factor is the same as the number of Feynman diagrams one can get if we consider 
different ways to contract the operators in the initial, intermediate and final states. For 
diagrams [T^)-d) Ns=8 ■ 2 ■ 4=. Here 8 comes from the eight ways to contract the nucleon 
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field operators in the initial and intermediate states; 2 is to account the interchange of two 
identical pions in final states; 4= is from the deuteron vertex. It should be noted that in this 
way the possibility to produce pions from both nucleon lines is included. The corresponding 
expressions of Fig. [T]d) can be obtained by changing C to "D in the above equations. 

The amplitudes of Figs. [1]:) and [T]i) can be obtained in a similar manner. The 
amplitudes of Fig. ^) are 



Af'^ = Ns- 1^^ iS^fe, ■ ,, ^' . ■ B^{1) , (A16) 



'"''^-^' ' -f, B',^ {V.,l) ■ -^ 

Itt^ r — pf — le 



■ , ^f . -BAip), (A17) 

p — Pi— IS 

j,c,pw _.r f l^dl p^ Mn f k^dk ■M^V^H/,^) 

BA{k) , (A18) 



k'^ — pf — ie 
j^cPW _ ,,_ ''^''' fp'dp A^^'^-^^^'fe,p) M^ ffdl ^p^,^^ ,^ Mr, 



ISI+In.SI - ^V5 ■ 2.^ J 2^2 (2M^)2 ■ p2_p2_,,- J ^^2 ^D [P, ^^ ^2 _ ^2 _ ,^ 



/ 



edk M}^k '(l^k) M 



N 



■BA{k). (A19) 



27r2 (2M^)2 P-pI- ie 

The corresponding amplitudes of Figs. [T]l) can be obtained by changing A to 5 in the above 
equations. 

2. Amplitudes of Fig. [l^)-h) 

For the diagrams of Fig. [I^)-h) there is no intermediate state interaction so that the 
amplitude can be written as: 

A^^ = A^,'^ + A^^,. (A20) 

Note that only the sum of the amplitudes of Figs. [Ife)-i) is independent of the choice of the 
pion field. That is why below we give only the amplitude for the sum of them. The partial 
wave amplitudes of Figs. [Ife)-i) without the ISI can be written as: 



Af'^ = Nl^^^{-4c, - 2c2 + J- - 2c3) ■ u{p) , (A21) 
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-Nl 



4M|m2 



P 

J n 



-4ci - 2c2 



9\ 



4M, 



N 



2C3) ■ W(p) , 



(A22) 



A 







^/2-^m| 



^4 

J 7r 



2^ 



u{l) 



dVtp dVti 



^,<M'i 



J n 



fdl 

2^ 



w{l) 



47r Att {p- If + ml 

dVLpdVLi l-3(p-f)2 
47r 47r (p'_r)2 + ^2 ^ 



(A23) 
(A24) 



.s+^+i,^^! ^ ^<,+ft+i7V2^iM^m2 /• /2c?/ |- dVtp dn^ 



^0 



"s 



m 



27r2 J 47r 47r ((p-- /)2 + ^2)2 



w(/) - ^ ) • (p - /")2 - 3wil) (pip ■ i) - C ' 



V2, 



Ag+h+i,'-^Di _ T^g+h+i'^gA^N^^-K [^ ^l f dVlp dVlf 



"S 



m 



2tt^ J An An 



MO 

V2 



(A25) 



{Ap'^-8p-l + 6l\l-pY -2f 



w{l) 



■ {f - /' + Ap- 1- 3p\p ■ If - 3l\p ■ If 



1 



(A26) 



((p-/)2 + m2)2 

The corresponding amplitudes with the ISI can be obtained through the following ex- 
pression: 



^•S'l.^^l r ^2 



m,PW Arm V^ f P dp M.]\f]\j 



PW 



27r2 



(2M;v) 



IN 

p2 — p| — ie 



^r'^'^'(p), (A27) 



where m can be 'e', 'f , or 'g+h+i'. The symmetry factors for these amplitudes are N^ = -^, 



Ni = ^ and Ar|+^+^ = ^. 



Appendix B: Observables 



In this appendix we present the expressions for the cross section near threshold. The 
amplitudes of pn — ?■ dTcn considered in this study are calculated at the threshold and, 
therefore, can be factored out of the phase space integration. The cross section for pn — )■ 
dn^TT^ is expressed in terms of the amplitudes given in Appendix A in the following way: 

"" = 4:2:2 ■ (^l-^/""!' + 31^/Sf) ■ ^EcMPCM 



(Bl) 
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where Eqm and pcM are the energy and momentum of the initial nucleons in the cm. frame. 
$ is the phase space factor defined as: 

*^/(2.)V(P-p..-p„-p.)p|^p|||-p||^ (B2) 

The prefactor ^^ is due to averaging over the initial spin states, the isospin wave function 
of the initial A^A^ states and the identity factor for the final pions, respectively. The factor 
of 3 in Eq. (IBip is from summing up the spin states. 
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